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ABSTRACT 

The relative extinction in the Galaxy computed with our new method (Choloniewski 
& Valentijn 1999, CV) is compared with three patterns: Schlegel, Finkbeiner & Davis 
(1998, SFD), Burstein & Heiles (1978, BH) and the cosecans law. It is shown that 
extinction of SFD is more reliable then that of BH since it stronger correlates with 
our new extinction. The smallest correlation coeffcient have been obtained for the 
cosecans law. Linear regression analysis show that SFD overestimate the extinction 
by a factor of 1.4. 

Our results clearly indicate that there is non-zero extinction at the Galactic South 
pole and that the extinction near the Galactic equator (|&| < 40°) is significantly larger 
in the Southern hemisphere than in the Northern. 

Key v^rords: dust, extinction - methods: statistical - Galaxy: general - galaxies: 
fundamental parameters 



1 INTRODUCTION 

We have introduced in our previous paper (Choloniewski 
& Valentijn 1999, hereafter CV) a new method for the de- 
termination of the extinction in our Galaxy. The method 
uses surface brightnesses of external galaxies in the B and R 
bands as listed in The Surface Photometry Catalogue of the 
ESO-Uppsala Galaxies (Lauberts & Valentijn, 1989, here- 
after ESO-LV) . The first draft of this method has been pub- 
lished in our earlier paper (Choloniewski & Valentijn 1991). 

The main purpose of the present paper is to compare 
our derived extinction values with recently published maps 
of extinction by Schlegel, Finkbeiner & Davis (1998), here- 
after SFD, and with the frequently used map of Burstein 
& Heiles (1978), hereafter BH (see also Burstein & Heiles 
1982). 



2 THE METHOD 

Our extinction determination (fully described in CV) em- 
ploys the surface brightnesses of external galaxies in the B 
and R bands: /is, fiR. Basically, our method produces the 
relative extinction compared to an overall mean - extinction 
with an unknown zero-point. The formula for the relative 
extinction (in B band) is simply a linear combination of /ib 
and fiR- 



Ab = 



1 — r s 



(1) 



In order to use equation (1) one has to know three parame- 
ters: r, s and c. 

The parameter r describes the ratio of extinction in 
the R and the B band (r = Ar/Ab) and is assumed to 
be constant. Its recent literature value is 0.61 (see CV for 
references) while we derived in CV two new estimates: 0.62 
and 0.64. As a reasonable compromise we adopt throughout 
this paper r = 0.62. 

The inverse of the parameter s describes the slope of the 
linear relation between surface brightnesses fis and fiR. The 
c parameter is introduced in order to maintain zero-point 
issues. Both parameters s and c depend on morphological 
type T, so they have been computed (using equations 7 and 
10 in CV) separately for every morphological type. 

In this paper we will consider several different subsam- 
ples. Extinction within every such subsample has been com- 
puted using the set of values s{T) and c(r) obtained from 
the same set of data. The s{T) and c(r) coefficients for the 
most important two subsamples used in ths paper are in 
Table 1. 

We use in this paper the extinction in B band as de- 
scribed in equation (1) and denote it as Ab{CV). 



3 THE SAMPLES 

For our analysis we use surface brightnesses at the radius of 
half total B light m B and R bands from ESO-LV. We exclude 
from the sample those galaxies which have morphological 
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Table 1. Coefficients s{T) and c(T) used for computing relative 
extinction in this paper according to equation 1 

sample "A" sample "B" 



T 


6 


i(T) 


c(T) 


s 


i(T) 


c(T) 


-5 


0, 


,826 


10.706 





.842 


10.233 


-4 


0, 


,955 


6.263 





,997 


4.199 


-3 


0, 


,915 


7.619 





,932 


6.927 


-2 


0, 


,893 


8.295 


0, 


,894 


8.306 


-1 


0, 


,785 


11.597 


0, 


,772 


12.060 





0, 


,768 


12.060 


0, 


,840 


10.129 


1 


0. 


,728 


12.972 





.737 


12.885 


2 


0. 


,691 


13.903 


0, 


.681 


14.166 


3 


0. 


,670 


14.223 





.668 


14.321 


4 


0. 


,652 


14.589 





.650 


14.663 


5 


0. 


,635 


14.979 





.637 


14.943 


6 


0, 


,641 


14.943 


0, 


.697 


13.627 


7 


0, 


,717 


13.163 


0, 


,714 


13.251 


8 


0, 


,731 


12.844 


0, 


.721 


13.181 


9 


0, 


,744 


12.656 


0, 


.784 


11.446 


10 


0. 


,795 


11.054 





.806 


10.675 



classifications suspected to be extinction dependent (marked 
in ESO-LV with Tfiag equal to 4). We also reject galaxies 
which have excessive (probably wrong) colours: hb — Hr 
smaller than zero or greater than 2.1. 

Since wc focus in this paper on the comparison of our 
extinction with SFD and BH data, we reject additionally 
those galaxies for which BH or SFD extinction values are 
not available. 

The largest complete subsamplc of galaxies of ESO-LV 
can be made by selecting galaxies which have visual apparent 
diameter Dorg greater or equal to 60 arcsoc. It contains, 
after applying the rejections described above, 7974 galaxies 
(sample "A"). 

As mentioned in CV, the completeness limit of ESO- 
LV galaxies is morphological type dependent and a univer- 
sal criterion, valid for every morphological type, is selecting 
galaxies with a visual diameter limit larger than 100 arcsec. 
We have used such more restricted sample for all the com- 
putations presented in CV. Wc use it also here as sample 
"B" (after rejections described above). 

Results described in Sections 4 and 5 show that more 
rigorousely defined sample "B" produces slightly more accu- 
rate extinction per galaxy than the sample "A" . However the 
sample "A" is three times larger than the sample "B" what, 
at least in part, compensate its slightly larger dispersion. 

Table 2 contains a summary of the definitions of sample 
"A" and sample "B". 

Since the ESO-LV galaxy catalogue covers the South- 
ern sky {5 < —17.5°) our analysis refers to this part of the 
hemisphere. 



4 THE ACCURACY 

Equation (1) represents an estimator of the foreground rel- 
ative extinction for an individual external galaxy. It is im- 
portant to know its uncertainty (standard error). In order 
to obtain this, we have divided the whole sky into squares 
with size A degrees and computed the average variance of 
the extinction inside these squares using the formula given 



Table 2. Definitions of the galaxy samples. 

General conditions: 

Ab{BH) - present 
Ab {SFD) - present 
Tflag 4 
0<jus -Mfl <2.1 

sample A sample B 

Dorg > 60 arcsec Dorg > 100 arcsec 
N=7974 N=2450 

by the so called one-way analysis of the variance theory (see 
e. g. Fish 1962). This average variance reflects the variations 
of the true extinction inside the squares with size A and the 
standard error of the extinction estimator expressed in equa- 
tion (1). So, when the size of the squares A tends to zero 
the average variance of the extinction should tend to the 
standard error of the extinction produced by our method. 

Figure 1 shows the average variance of extinction as a 
function of A. The minimum value for A which we applied 
was one degree. For smaller values of A a too large fraction 
of studied squares contain only one galaxy (such squares 
can not be taken into account in the one-way analysis of 
variance) . 

As we expect, the average variance has minimum for 
the smallest applied value of A (one degree). For sample 
"A" the standard error of our extinction in Ab is 0.43 mag- 
nitude, while for sample "B" this is 0.40 magnitude, which 
corresponds to a standard error in E{B — V) of approxi- 
mately 0.10 magnitude (for Ab/E{B - V) = 4.3). Sample 
"B" produces relative extinction with slightly higher accu- 
racy than sample "A". 

BH used for the calibration of their extinction map B- 
V colours of 131 globular clusters and RR Lyrae stars. SFD 
used for the calibration of their extinction map B-R colours 
of 106 brightest cluster ellipticals and B-V colours of 389 
elliptical galaxies with measured Mg2 index (505 objects in 
total). 

Both BH and SFD report that their calibrators show a 
residual scatter in B-V, with respect to calibration regression 
line, of approximately 0.03 magnitude. 

Our extinction estimator, when applied to photographi- 
cally measured surface brightnesses of galaxies in two bands 
as listed in the ESO-LV catalogue, has three times larger 
standard error than the calibrators used by BH and SFD. 
This is an important disadvantage (at least as long as we 
apply it to ESO-LV photometrical data). But there is one 
important advantage of our extinction estimator - it can 
be applied to many more objects since wo can apply the 
method, at present, to 7974 galaxies (sample "A") from the 
ESO-LV catalogue. 



5 THE CORRELATION 

We have computed the Pearson, Spearman and Kendall cor- 
relation coefficients (see Press, Teukolsky, Vetterling & Flan- 
nery 1992 for definitions and software) between our extinc- 
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Table 3. Pearson, Spearman and Kendall correlation coefficients 
between Ag^CV) extinction estimate and three other extinction 
estimates: SFD, BH and csc(b). 

Pearson Spearman Kendall 
Sample "A" 



SFD 0.276 

BH 0.247 

csc(b) 0.202 

Sample "B" 

SFD 0.327 

BH 0.293 

csc{b) 0.238 



0.234 
0.219 
0.188 



0.304 
0.282 
0.252 



0.158 
0.148 
0.127 



0.208 
0.192 
0.170 



tion and extinction given by SFD and BH and for the cose- 
cans law: 



Ab = Ao CSC 



(2) 



where b denotes Galactic latitude. The computation have 
been performed for sample "A" and sample "B" (see Ta^ 
ble 3). The parameters s{T) and c(T) have been computed 
separately for every sample (sec Table 1). 

All three correlation coefficients for both samples are 
the largest for SFD extinction and the smallest for the cose- 
cans law. The extinction of BH is always between these two 
extreme results. 

The coefficients are generally higher for sample "B" 
than for sample "A" what suggests that sample "B" pro- 
duces more accurate extinction than sample "A". 

Since the sample "B" is a subsample of the sample "A" 
the correlation coefficients for "A" and "B" are not statisti- 
cally independent. In order to produce a set of statistically 
independent correlations we divide the sample "A" into six 
subsamples: 

(i) 60 arcsec < Dorg < 80 arcsec (N=3876) 

(ii) %Q arcsec < Dorg < WO arcsec (N=1648) 

(iii) 100 arcsec < Dorg < 120 arcsec (N=718) 

(iv) 120 arcsec < Dorg < 140 arcsec (N=647) 

(v) 140 arcsec < Dorg < 160 arcsec (N=298) 

(vi) 160 arcsec < Dorg (N=787) , 

and compute the parameters s{T) and c(T) separately for 
every subsample. As a result of this procedure we have six 
statistically independent correlation coefficients - see Figs 2, 
3 and 4 for results. As before, the correlations are the largest 
for SFD, smaller for BH and the smallest for cosecans law. 



6 CORRECTION OF THE BH'S EXTINCTION 
ZERO POINT 

The formulae of BH produces for some regions of the sky 
extinction less than zero (for the samples analyzed in this 
paper the minimum value of Ab{BH) extinction is -0.12 
magnitude). In spite of BH's instructions to set these val- 
ues to zero we have actually used these negative values and 
found that Ab(CV) is for them significantly smaller than 
for Ab{BH) « (see upper panel of Figs 7 and 8). This 
means that the BH extinction (in the B band) was underes- 
timated by 0.12 magnitude - just the absolute value of the 



minirrmm value of Ab{BH) extinction. This is in approxi- 
mate agreement with SFD who discovered a similar offset 
of 0.09 magnitude. We use in the regression analysis pre- 
sented in the next Section the corrected BH's extinction: 
Ab{BH)c = Ab{BH) + 0.12 instead of Ab{BH) itself. 



7 THE LINEAR REGRESSION 

In the ideal case there would be linear dependence with slope 
equal to one between our relative extinction estimate (CV) 
and the extinction of SFD and BH (corrected). 

Since our extinction values are relative, with an arbi- 
trary zero-point, the constant term of this linear dependence 
should not be equal to zero. The constant, multiplied by —1, 
should be added to our relative extinction to transform it to 
the absolute extinction. 

We have fitted the straight lines (using the least squaxes 
method) taking as independent variables extinction of SFD 
and BH and as dependent variable our estimate of relative 
extinction (CV). We have found the following regression co- 
efficients for sample "A": 

A B (CV) =0.662 (±0.028) As (SFD) - 0. 189 (±0.010) (3) 
Ab{CV) = 0.555 {±0.026)Ab{BH)c - 0.176 (±0.010) (4) 
and for sample " B" : 

ABiCV) = 0.741 (±0.043) As (SFD) - 0.225 (±0.016) (5) 

Ab{CV) = 0.608 {±0M0)Ab{BH)c - 0.206 (±0.017) (6) 

where in brackets la errors are given. The differences be- 
tween slopes and constant terms computed for sample "A" 
and "B" are only marginally larger than the combined er- 
rors. In the forthcoming we use their averages. 

A graphical presentations of the regression lines are 
shown in Figs 5-8. In order to investigate whether the pos- 
tulate about linear dependence is valid, the regression lines 
are shown together with row data (lower panels) and with 
averages of Ab{CV) computed for 0.05 magnitude bins of 
ABiSFD) and Ab{BH) (upper panels). 

The slopes of the regression lines are definitely less than 
one: for SFD the slope is circa 0.7 while for BH it is circa 0.6. 
The constant terms for both SFD and BH are approximately 
the same: -0.20 magnitude. This defines the zero point of 
our results: adding 0.20 magnitude to our relative extinction 
transforms it to the absolute extinction. 

When we introduce our absolute extinction into equa- 
tions 3-6, the constant terms in all 4 equations gets close 
to zero, now setting the relation between our (absolute) ex- 
tinction and SFD's and BH's extinction. Thus we conclude: 
SFD overestimate extinction by a factor of 0.7~^ w 1.4 while 



BH by a factor of 0.6 



1.7. 



Up to now five papers reported similar results, namely: 
that SFD extinction is about 1.4 times too large. Stanek 
(1998b) obtained, using colours of low galactic latitude glob- 
ular clusters, that the overestimation factor is 1.35 (see 
also Stanek 1998a). Arce & Goodman (1999a, 1999b) an- 
alyzed extinction in the Taurus region using four indepen- 
dent methods and found that the factor is 1.3 - 1.5. Gonza- 
lez, Fruchter & Dirsch (1999) obtain extinction in the small 
field around GRB 970228. They found using two methods 
Av = 0.55 while Av{SFD) = 0.78. The ratio of these two 
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numbers give again a factor 1.4. Our finding that SFD over- 
estimate extinction is also supported by Ivans et al. (1999) 
and von Braun & Mateo (2001). 

Some kind of doubt about the correctness of the cahbra- 
tion is also in original SFD paper where wc can find as the 
comment to their Fig. 6 the statement that "A slight trend 
in the residuals is evident for both BH and DIRBE/IRAS 
corrections, in the sense that the highest reddening values 
appear to be overestimated" . 

Inspecting Figs 5-8 demonstrate a general linear depen- 
dence between our relative extinction estimate and the ex- 
tinction of SFD and BH. However, some minor, but statis- 
tically significant, deviations arc clearly present, especially 
for Ab{SFD) and Ab{BH) less than 0.5 magnitude. 



8 GALACTIC LATITUDE DEPENDENCE 

Unfortunately, wc can not use our method of extinction de- 
termination, as applied to the ESO-LV galaxy catalogue, to 
produce a new high resolution extinction map and to com- 
pare it with SFD and BH maps. This is because our data 
are too sparse (less than one galaxy per square degree) and 
their accuracy is too low {a(E{B — V)) ~ 0.10 per galaxy) 
mainly because of instrinsic scatter of the surface brightness 
of galaxies. 

But the quality of our data is good enough to evaluate 
their galactic latitude dependence. We show this dependence 
for our absolute extinction (equal to the relative extinction, 
as defined in equation 1, plus 0.20 magnitude) together with 
SFD, corrected BH and the cosecans low - see Fig. 9. All four 
solutions have been averaged inside 5 degrees galactic lati- 
tude bins in the points where ESO-LV galaxies are (sample 
"A'"). Wc use Ao = 0.10 for the cosecans law (equation 3) 
to match as close as possible to our solution. 

Fig. 9 shows that our absolute extinction is always 
greater than zero confirming the value of the zero-point com- 
puted in Section 7. 

As one can see in Fig. 9 our solution as well as SFD and 
BH mimic quite well the classic cosecans law and indicate 
considerable extinction near the galactic pole. The differ- 
ences between our results and SFD and BH are the largest 
near the galactic equator and at the northern hemisphere. 

Since we find that both SFD and the (corrected) BH 
extinction standards overestimate Galactic extinction by a 
factor 1.4 and 1.7 respectively, we show the same data as 
in Fig. 9 but with SFD and corrected BH extinction divided 
by these factors (or, equivalently, multiplied by 0.7 and 0.6 
respectively) - see Fig. 10. The better agreement of our so- 
lution with the rescaled SFD and BH data is evident, es- 
pecially near the Galactic equator (fe ~ 0), confirming the 
need for rescaling SFD and BH extinction. But even for the 
rescaled data statistically significant differences between our 
results and SFD and BH, although smaller, still exist. Es- 
pecially our data exhibit a South - North asymmetry near 
the galactic equator (|6| < 40°) with more extinction in the 
Southern Galactic hemisphere which is not visible in SFD 
and BH data. 



9 DISCUSSION 

Our extinction estimator (introduced in CV) definitely indi- 
cates that the new extinction map of SFD is more realiable 
than the old one of BH. The historical cosecans law is in this 
competition on the last place. This result has been obtained 
using six separate, statistically independent, sets of data. 

The superiority of the SFD extinction map over the 
BH map have been demonstrated by computing correlation 
coefficients with our extinction results. However, the cor- 
relation coefficients are not sensitive to the amplitude of 
variation of the input data nor its zero-point. In order allow 
for further, more specific, comparisons we have performed 
a linear regression analysis between our relative extinction 
and SFD's and BH's extinction. This analysis provides us 
the zero-point of our extinction (equal to 0.20 magnitude) 
which allows us to transform our relative extinction into ab- 
solute extinction. We have showed that, in comparison to 
our absolute extinction, SFD overestimate extinction by a 
factor of 1.4. This is in agreement with five other authors 
and considerably changes our view on the amplitude of the 
Galactic extinction. 

Superiority of SFD over BH reddening map, as we re- 
port in this paper, does not mean that the first is an ideal 
result - we have discovered some significant differences be- 
tween our extinction and SFD when analyzing their Galactic 
latitude dependence even after correcting SFD extinction by 
dividing it by 1.4 factor. One possible source of these dif- 
ferences is that SFD assume that the extinction is strictly 
proportional to the dust column density what need not be 
true - the extinction to dust ratio can vary across the sky 
and may also depend on the wavelength. 

Our results are important for creating reliable extinc- 
tion standards and demonstrate the correctness and use- 
fulness of our extinction estimator. They can be also re- 
garded as a stimulus for applying our extinction estimator 
to other, larger and more accurate galaxy catalogues, partic- 
ularly thoose expected from new wide field imaging surveys, 

Future applications of our method to other databases 
may produce maps with a courser resolution. At present the 
method applied to ESO-LV data is unable to generate the 
extinction map with both sufficient resolution and accuracy. 
Our present method should be rather taken as new extinc- 
tion calibrator. But using the method to larger and more 
precise galaxy catalogues may result into stand-alone high 
resolution extinction maps. 
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Figure 1. The average variance of extinction Ab{CV) inside A degrees squares on the sky as a function of the size (A) of these squares. 
The minimum variance (at one degree: iop(A) = 0) represent the standard deviation of the estimator of relative extinction used in this 
paper. 
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Figure 2. Pearson correlation coefficient between our Ab{CV) extinction and the extinction of SFD and BH and the cosecans law 
computed for six subsamples defined using visual diameter Dorg (see Section 5 for details) 
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Figure 3. The same as Fig. 2 but for Spearman correlation coefficient. 



Extinction in the Galaxy from surface brightnesses of ESO-LV galaxies: testing 'standard' extinction maps 9 



.5 



1 r 



.3 



Kendall 



.1 



: ' / \ 



SFD 

BH 

csc(b) 







50 



100 



150 



D 



org 



J L 



200 



Figure 4. The same as Fig. 3 but for Kendall correlation coefficient. 
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Figure 5. The dependence between the extinction of SFD Ab(SFD) and our relative extinction Ab(CV). Upper panel shows the 

averaged data inside 0.05 magnitude bins. Error bars represent standard deviation (lo"). Lower panel shows raw data. The data axe taken 
from sample "A" (see text). Dotted straight lines visible on both panels represent the least squares fit. 
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Figure 6. The same as Fig. 5 but for sample "B". 
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Figure 7. The same as Fig. 5 but for extinction of BH Ab{BH). 
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Figure 8. The same as Fig. 7 but for sample "B". 
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Figure 9. Galactic latitude dependence of our absolute extinction (represented by the solid broken line with Icr error bars) compared 

with the extinction according to SFD, BH and the cosecans law. The BH extinction have been corrected by adding 0.12 magnitude 
constant. Note north-south assymetry of our extinction near the Galactic equator (|b| < 40°) and non-zero extinction near the Galactic 
south pole (6 » —90). 
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Figure 10. The same as Fig. 9 but with SFD extinction multiplied by the factor 0.7 and for corrected BH extinction multiplied by 0.6. 



